home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / dlaqsb.z / dlaqsb
Text File  |  1998-10-30  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAQQQQSSSSBBBB((((3333FFFF))))                                                          DDDDLLLLAAAAQQQQSSSSBBBB((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAQSB - equilibrate a symmetric band matrix A using the scaling factors
  10.      in the vector S
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAQSB( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED, UPLO
  16.  
  17.          INTEGER        KD, LDAB, N
  18.  
  19.          DOUBLE         PRECISION AMAX, SCOND
  20.  
  21.          DOUBLE         PRECISION AB( LDAB, * ), S( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      DLAQSB equilibrates a symmetric band matrix A using the scaling factors
  25.      in the vector S.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              Specifies whether the upper or lower triangular part of the
  31.              symmetric matrix A is stored.  = 'U':  Upper triangular
  32.              = 'L':  Lower triangular
  33.  
  34.      N       (input) INTEGER
  35.              The order of the matrix A.  N >= 0.
  36.  
  37.      KD      (input) INTEGER
  38.              The number of super-diagonals of the matrix A if UPLO = 'U', or
  39.              the number of sub-diagonals if UPLO = 'L'.  KD >= 0.
  40.  
  41.      AB      (input/output) DOUBLE PRECISION array, dimension (LDAB,N)
  42.              On entry, the upper or lower triangle of the symmetric band
  43.              matrix A, stored in the first KD+1 rows of the array.  The j-th
  44.              column of A is stored in the j-th column of the array AB as
  45.              follows:  if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-
  46.              kd)<=i<=j; if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for
  47.              j<=i<=min(n,j+kd).
  48.  
  49.              On exit, if INFO = 0, the triangular factor U or L from the
  50.              Cholesky factorization A = U'*U or A = L*L' of the band matrix A,
  51.              in the same storage format as A.
  52.  
  53.      LDAB    (input) INTEGER
  54.              The leading dimension of the array AB.  LDAB >= KD+1.
  55.  
  56.      S       (output) DOUBLE PRECISION array, dimension (N)
  57.              The scale factors for A.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAQQQQSSSSBBBB((((3333FFFF))))                                                          DDDDLLLLAAAAQQQQSSSSBBBB((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      SCOND   (input) DOUBLE PRECISION
  75.              Ratio of the smallest S(i) to the largest S(i).
  76.  
  77.      AMAX    (input) DOUBLE PRECISION
  78.              Absolute value of largest matrix entry.
  79.  
  80.      EQUED   (output) CHARACTER*1
  81.              Specifies whether or not equilibration was done.  = 'N':  No
  82.              equilibration.
  83.              = 'Y':  Equilibration was done, i.e., A has been replaced by
  84.              diag(S) * A * diag(S).
  85.  
  86. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  87.      THRESH is a threshold value used to decide if scaling should be done
  88.      based on the ratio of the scaling factors.  If SCOND < THRESH, scaling is
  89.      done.
  90.  
  91.      LARGE and SMALL are threshold values used to decide if scaling should be
  92.      done based on the absolute size of the largest matrix element.  If AMAX >
  93.      LARGE or AMAX < SMALL, scaling is done.
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.